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A non-conforming domain decomposition approximation for the 
Helmholtz screen problem with hypersingular operator * 
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Abstract 

We present and analyze a non-conforming domain decomposition approximation for a 
hypersingular operator governed by the Helmholtz equation in three dimensions. This oper¬ 
ator appears when considering the corresponding Neumann problem in unbounded domains 
exterior to open surfaces. We consider small wave numbers and low-order approximations 
with Nitsche coupling across interfaces. Under appropriate assumptions on mapping prop¬ 
erties of the weakly singular and hypersingular operators with Helmholtz kernel, we prove 
that this method converges almost quasi-optimally. Numerical experiments confirm our error 
estimate. 

Key words: Helmholtz problem, hypersingular operator, boundary element method, domain 
decomposition, Nitsche method. 
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1 Introduction 

In recent years we have started to develop non-conforming boundary elements, in the sense 
that approximations to boundary integral equations with hypersingular operators can be dis¬ 
continuous. Approaches consider both element-wise discontinuous methods [mill! and domain 
decomposition techniques, mortar coupling in m and Nitsche coupling in However, all 
results are restricted to the simple model problem of the Laplacian. 

In this paper we extend the Nitsche domain decomposition method from |1| to the hypersin¬ 
gular operator Wk stemming from the Helmholtz problem with small wave number k. Traditional 
variational analysis of this operator is based on the theory of Fredholm operators since, for small 
wave numbers, Wk can be handled as a compact perturbation of the elliptic operator Wq which 
corresponds to the Laplace case. This approach is not applicable to our non-conforming discrete 
setting. The energy space of Wq, e.g. defined on an open surface T, is a trace space of H^{Q) 
(with n := M^\f) and thus of order 1/2. In such a space there is no well-defined trace operator. 
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On the other hand, the analysis of discontinuous approximations requires the consideration of 
jumps and thus, traces. Because of this conflict, numerical analysis of discontinuous approx¬ 
imations of hypersingular integral equations has been carried out exclusively on the discrete 
level where traces are defined as restrictions. In this way arguments from variational settings 
can be avoided. Now, standard numerical analysis of Fredholm operators is based on compact¬ 
ness arguments which, by nature, are connected with non-discrete variational settings, cf., e.g., 
p3l fT6| fT8l C] where the analysis of boundary elements is based on Garding’s inequality. In this 
paper, we present an analysis of the Helmholtz case which reconciles both seemingly conflicting 
approaches, the restriction to discrete spaces and appropriate extension to consider Fredholm 
operators. This latter extension is done by providing discrete variants of a Garding’s inequality. 
Nevertheless, our main result will be based on three assumptions on the weakly singular and 
hypersingular operators whose verihcation goes beyond the scope of this paper. 

Our analysis also uses a compactness argument. Corresponding estimates generate unknown 
constants which depend in most cases on the geometry and possibly other data; in our case they 
depend on the order of Sobolev norms. For this reason, final estimates are based on Sobolev 
regularities s > 1/2. Limits of s tending to 1/2 cannot be considered since the dependence of 
the constants on s is unknown. This is different in the Laplace case where estimates involving 
natural norms of order 1/2 can be established by limits. In this way quasi-optimal error estimates 
with poly-logarithmic perturbations appear, cf., e.g., [H Theorem 3.1]. In the Helmholtz case 
considered here, estimates are less specific by assuming that Sobolev orders in upper bounds are 
strictly larger than 1/2. 

Let us note a further complication of discontinuous boundary elements. Discontinuous (DG) 
finite elements are usually analyzed considering specihc DG-type norms, comprising broken 
semi-norms and scaled jump terms. They are tuned to harmonize with DG-bilinear forms and 
have also been considered in the boundary element settings studied in [Enani!. However, in 
a boundary integral operator approach one has to consider a post-processing step consisting in 
evaluating the underlying representation, e.g., Uh{x) = opp(iL, u/i) for x G H {uh denoting the 
boundary element approximation, and opp the integral operator with kernel K used for repre¬ 
senting the solution U = opp(iL,n) to the original boundary value problem). One establishes 
convergence orders for the point-wise evaluation by applying duality estimates to the integral 
operator, 

\U{x) - Uh{x)\ = |opp(iL,u - Uh)\ < ||u - Uh\U\\K{- - x)||*/. (1) 

Here, || • ||* and || • ||*/ denote, respectively, the norm considered to bound the boundary element 
error u — Uh and its dual norm, and one uses that the kernel K{y — x) is smooth for x / y (x G H, 
y G F). It is not straightforward to analyze such a duality estimate for a DG-norm. In this paper 
we provide an error estimate for a standard Sobolev norm (it is a broken Ll^/^-norm) so that its 
dual norm is known and can be used to make the error estimate 0 explicit by specifying both 
norms. 

The remainder of this paper is organized as follows. In the next section we briefly recall 
some Sobolev norms and present the model problem. We also formulate two assumptions on 
which our subsequent analysis is based. In Section we present the non-conforming domain 
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decomposition setting and formulate the main result (Theorem]^, a Cea-type estimate. The 
following Corollary establishes the convergence order of the method. A proof of Theorem 
is given at the end of Section after collecting a number of preliminary results, including 
consistency of the discrete method (Lemm a [^ , boundedness of the sesquilinear form in broken 
Sobolev spaces of order s > 1/2 (Lemma 10), discrete Garding’s inequalities (Lemma 11 and 
Corollary 12), and a lower-order error estimate based on the Aubin-Nitsche trick (Lemma 13). 
Some numerical experiments that confirm our estimates are reported in Section 

Throughout the article, we will use the symbols and in the usual sense. In short 
ah{v) < bh{v) when there exists a constant C > 0 independent of v and the mesh size h, such 
that ahiv) < Cbh{v). Also, ahiv) ~ bhiv) means that ativ) < bh{v) and ah{v) > bh{v). 


2 Sobolev spaces and model problem 

For 0 C M"' and 0 < s < 1 we define 


2 



+ \u 


2 

Hop) 


with semi-norm 


^Hop) 



u{x) - u{y)\^ 
'~\x - y|25+n- 


1/2 


( 2 ) 


For a Lipschitz domain 11 and 0 < s < 1, the space H^{Q) is defined as the completion of (7^(11) 
under the norm 

II II ft |2 f j 

Mfiop) ■- dist(x,5L!)2- • 

For s G (0,1/2), || • and || • are equivalent norms whereas for s G (1/2,1) there 

holds i7®(0) = Hq{Q), the latter space being the completion of (7“ (11) with norm in Lf^(ll). 
For s > 0 the spaces and if“*(11) are the dual spaces (with L‘^{Q) as pivot space) of 

17*(0) and 77*(H), respectively. For more details on Sobolev spaces we refer to |17l llOj . 

In the following, let F be a piecewise plane Lipschitz surface. For simplicity we assume that 
F is open with polygonal boundary clF. Sobolev spaces on faces of F are defined as previously, 
identifying faces with sub-domains of i.e., n = 2 in (§. For a closed surface F being the 
boundary of 11 and containing F, 71*(f) is the trace of {s > 1/2) and H^(T) is the 

space of functions from 77* (F) with support on F. Dualities with spaces of negative order are 
defined as previously. Furthermore, throughout the paper, we use the same notation for Sobolev 
spaces of vector-valued functions, taking respective norms component-wise. 

Our model problem is: For given wave number k > 0 and sufficiently smooth function f find 
u G H^f‘^{T) such that 


Wku{x) := 


1 d 
Att dn^ 


u{y) 


Q Qik\x-y\ 
diiy \x - y\ 


dSy = f{x), 


X G F. 


Here, n is a normal unit vector on F pointing to one side. 


( 3 ) 
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Remark 1. The assumption that F is an open surface implies that (|^ has a unique solution 
for any k with Im{k) > 0, cf. For closed surfaces, different boundary integral equations 

like the Burton-Miller formulation are in order, see m- 

Remark 2. (i) In the case of the Laplacian, i.e., k = ^, it is well known that the solution of 
Q with appropriate (and sufficiently smooth) right-hand side f (so that it relates to a Neu¬ 
mann Laplace problem) satisfies u G iF^(r) for any r < 1, see \^. In the Helmholtz case 
(k > d) Stephan used the theory of pseudo-differential operators to show that on open surfaces 
with smooth boundary curve, and f G iF^(r), u has a square-root edge singularity and that 
u G W^iV) for any r < 1. We do not know of a specific analysis on open or closed polyhedral 
surfaces. 

(ii) A direct formulation of the Helmholtz problem in \ F with Neumann boundary condi¬ 
tion satisfies the Sommerfeld radiation condition, i.e., it only considers outgoing waves. The 
boundary integral equation with hypersingular operator Wk for the wave number k reflects this 
behavior. Changing the sign of k turns the problem into the non-physical one of incoming waves. 
In our analysis we will need the adjoint operator ofWk. It can be immediately seen that this 
is W-k, when considering the lf{T)-sesquilinearform. Therefore, mapping properties of W-k 
can be proved analogously to the ones of Wk by replacing the Sommerfeld radiation condition of 
outgoing waves by the one representing incoming waves, cf, e.g., and see also pM Remark 
3.9.6]. However, for the particular case of an open polyhedral surface the literature is scarce, as 
most specific results concern the Laplacian. 

Considering the two previous remarks, we are making the following assumptions. 

Assumption 1. There exists r^eg G (1/2,1) such that the solution u of i) satisfies u G 

#Ueg(F), 

Assumption 2. There exists r\so £ (1/2,1) such that, fork > 0, the operator W-k : R'’‘=°(F) —> 
iF'’i=°-^(F) is an isomorphism. 

A variational formulation of ^ is; Find u G R^/^(F) such that 

{Wku, v)r = if, v)r Vn G H^/\r). (4) 

Here, (•, •)r denotes the duality pairing between R~^/^(F) and R^/^(F). Throughout, this generic 
notation will be used for the L^-inner product and other dualities, and the domain is indicated 
by the index. 

A standard boundary element method for the approximate solution of Q is to select a 
piecewise polynomial subspace Hh C H^^‘^{T) and to define an approximant Uh G Hh by 

{WkUh, v)r = if, v)r Vu G Hh. 

3 Domain decomposition with Nitsche coupling 

In this section, we introduce the Nitsche-based boundary element method for the approximate 
solution of problem Q, and present the main result. Theorem]^ 
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3.1 Preliminaries 

We consider a decomposition of F, 


T ■■= {Fj-; j = 1,..., J}, 


where we assume that elements of T are plane polygonal surfaces. Throughout the paper, we 
will use the notation Vj for the restriction of a function u to a sub-surface Tj (also called sub- 
domain). The decomposition of F induces product Sobolev spaces of complex-valued functions, 
e.g., 

H^{T) := 

with corresponding broken semi-norm |•|Hs( 7 '), using on each sub-domain the Sobolev-Slobodeckij 
semi-norm previously defined. This notation with decomposition T will be used generically, i.e., 
also for the piecewise L^-sesquilinear form 

{v,w)r ■■= 

j 


and its extension by duality to H^{T) x 

We also make use of the surface differential operators curl and curl. On a subset of x {0} 
they amount to curl(^ := O) and curlc^ := dxi^P 2 — dx^^fi for sufficiently smooth 

scalar and vector functions tp and cp = {ipi,ip 2 ,P 3 )'^, respectively. For a definition and analysis 
on Lipschitz surfaces we refer to [2j. The restrictions of these operators to a face Fj will be 
denoted by curlj and curlj. Corresponding to the decomposition 7~ we also define the broken or 
piecewise operators cufIt- and curlT-, e.g., (cufIt-:= cufIj pj (j = 1,..., J) and similarly 
the other operator. 

Let 7 denote the skeleton of T, including dT. The jump [•] of functions across 7 is defined so 
that it is compatible with a tangential direction on 7 , appearing when integrating by parts the 
surface differential operators. More precisely, for a scalar function v (sufficiently 7~-piecewise 
smooth) and a tangential vector field p (sufficiently smooth so that its trace on 7 is well defined) 
we define tangential components t((^) and jumps [u] being compatible with the integration-by¬ 
parts formula 

H )7 = (cuiir p, v)r - (cuFlr u, p)t- (5) 

We select a unique tangential direction on 7 \ (9F (this fixes the directions of the jumps), and 
on ST so that [u]|ar is the trace of u on 5F. 

Now, for s G [1/2,1] and > 0, we introduce the norm 

\\'^\\ h ={ t ) ■= (l^lH‘>(r) + 

For s > 1/2, this is a norm in H^{T) and in the case s = 1/2, this norm will be used only for 
discrete functions whose jumps across 7 are well defined as elements of ^^( 7 ). 
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We end this section with recalling a relation that connects the hypersingular operator Wk 
with the single layer operator I4 defined by 

H n 2fc I IE y\ 

14 (^(x) := — / ip{y)- - rdSy, x G T. 

dvr Jy \x-y\ 

When applied component-wise to vector-valued functions we use the bold face symbol T4- The 
operators Wk and I4 satisfy the relation 

(Wfctt, u)r = (Vfc curlu, curlu)r — fe^CVfcntt, nu)r Vu, u G (6) 

see naEi]. As in previous publications on the Laplacian, this formula will give rise to our 
non-conforming discrete formulation of the hypersingular operator. 

3.2 Discrete method and main result 

On every sub-domain Fj we consider regular, quasi-uniform meshes 7j, j = 1,..., J, of shape- 
regular elements (quadrilaterals or triangles), = iJKeTjK. The maximum, respectively mini¬ 
mum, diameter of the elements of Tj is denoted by hj, respectively hj. We also define 

h := max{/ii,..., /ij}, h := min{/i^,..., hj}. 

Throughout this paper we assume that 0<h<h<C< 00 . Indeed, our main result assumes 
globally quasi-uniform meshes (h — h). But since some technical results hold for more general 
meshes we use the notation of h. We introduce discrete spaces on sub-domains consisting of 
piecewise (bi)linear functions: 

Xkj := {u G C^Tj); v\k is a polynomial of degree one ViF G 7)}, j = 1,..., J. 

Our global approximation space then is 


Xh TljXhj- 

We identify both product spaces H^{T) and X^ with their direct sums, e.g., X^, = X^p © • • • © 
Xh^j so as to consider their elements as scalar functions. Doing so, we note that Xh {Zi 
due to the possible discontinuity and non-vanishing trace on dT of its elements. Using the 
discrete space Xh for the approximation of ^ requires a different sesquilinear form that is well 
defined for such functions and that controls their jumps. 

For given > 0 and r G M, we dehne the following sesquilinear form on Xh x Xh- 

Ar(v, w) := {Vr cufIt- V, curly- w)-] — r^(Wnu, ntc)r 

+ {Triv), + {[v],T_riw))^ + I^{[v], [u;])^ 

with operator being given by (cf. @) 

Tr{v) := t(W curly- V G H^{T), s > 1/2, r G M. (7) 
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The Nitsche-based non-conforming domain decomposition method associated to problem Q 
then reads as: Find Uh G such that 

Ak{uh,v) = {f,v)r 'iv € Xh- (8) 

The analysis of this scheme will be based on a third assumption which is quite natural but 
whose proof we have not found in the literature for our precise situation. It is well known that 
14 and Wk are Fredholm operators of index zero. This follows from the fact that they are, 
respectively, compact perturbations of the positive definite operators Vq and Wq as mappings of 
their energy spaces to the dual spaces. For a closed smooth surface this follows from the theory 
of pseudo-differential operators and has been extended by Stephan [23] to open surfaces. There, 
it is shown that 

G (0,1/2] : Vk - Vq : ^ (9) 

for a smooth open surface. Our assumption is that this holds for our open, piecewise plane 
Lipschitz surface. 

Assumption 3. There holds (§. 

The main result of this paper is: 

Theorem 3. Let Assumptions 1,2,3 hold true and assume that the meshes defining Xh are 
globally quasi-uniform, i.e., h h. Given e > 0 choose v ~ h~'". Then the discrete scheme Q 
is uniquely solvable for h small enough. Furthermore, selecting e' > e and s G (l/2,rreg], there 
exists ho > 0 such that there holds the almost quasi-optimal error estimate 

||ri - Uh\\Hi/ 2 (r) ^ inf jjtt - v\\hs(t) < ho. 

Here, u and Uh are the solutions of Q and Q, respectively. 

A proof of this result will be given at the end of Section [^ We also obtain the following a 
priori error estimate. 

Corollary 4. Let Assumptions 1,2,3 hold true and assume that the meshes defining Xh are 
globally quasi-uniform. Given e' > e > 0 choose v ~ Then there exists ho > 0 such that 

there holds 

||tt - V/l < ho. 

Proof. We combine the error estimate by Theorem with standard approximation properties. 
The assertion follows by selecting s = 1/2with e" > 0 and renaming 2e''-|-e' as a new eh □ 


4 Technical details and proof of the main theorem 


We start with collecting some preliminary technical results in the following subsection. Then, 
we prove essential ingredients of the proof of Theorem which is given at 


in Subsection 4.2 


the end of this section. 
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4.1 Preliminary results 

We will make use of the continuity (see 0): 


Vr : 0 < s < 1, re 


( 10 ) 


Proofs for the statements of the following lemma can be found in |12[ Lemma 5] and [9l Lemma 
4.3], 

Lemma 5. Let R C be a Lipsehitz domain with boundary dR. 

(i) There holds 

1 „ „ 

( 11 ) 


^ 1/2-|s| 


for any s e (—1/2,1/2) and any v e H^{R). 
(a) There holds 


’Wl^idR) 


< 


1 




( 12 ) 


(13) 


for any s e (1/2,1] and any v e H^{R). 

Lemma 6. For r e M there holds 

\\Trv\\L^^) < (s - l/2)"^/^|w|ffo(r) Vr e H^{T), 1/2 < s < 1, 
with hidden constant depending on r. 

Proof. For the case r = 0, this estimate has been shown in Lemma 4.2], Using the continuity 
(10) for wave number r ^ 0, the same estimates apply. □ 

(14) 

(15) 


Lemma 7. There holds 


curlru||^-i/ 2 ( 7 -) < (s - 1/2) ^\v\Hs{r) Vr e H^{T), l/2<s< 1 , 
II curlr r||j:^-i/ 2 ( 7 -) - \'>^\m/2{T) Vr e 


Proof. The first estimate can be proved by using the equivalence ( [IT| ) of H^{Tj) and H^fTj)- 
norms for s G ( — 1/2,0], the continuity of curlpj : H^{Tj) —)• Lf^^^Fj) {s G (1/2,1]) and a 
quotient-space argument (cf. [I3]). For details see [U (4.15)]. Bound (15) follows by applying 
face-wise the corresponding estimate from la Lemma 4.1]. □ 


Lemma 8. For n > 0 and s > 1/2 there holds 

ll^llL 2 (r) < \\v\\h^{T) Vr G H^{T). (16) 

Proof. The proof is a slight variation of the proof of the Poincare-Friedrichs inequality by a 
compactness argument. In this case, we use the compactness of the embedding of H^{/T) in 
L^(F) and the boundedness of the functional ([•], [•])y^ on H^i/T). Furthermore, the kernel 
of j • \h‘’(T) consists of T-piecewise constant functions which are eliminated by the functional 

([•], [•])-y' (note that the jump [•] reduces to the trace operator on (9F). □ 







4.2 Consistency, boundedness, discrete ellipticity, and Aubin-Nitsche trick 


In this section we show four essential ingredients for the proof of Theorem These are the con¬ 
sistency of the non-conforming discrete scheme (Lemma |^, the boundedness of the sesquilinear 
form (Lemma 10), its discrete ellipticity in the form of Carding inequalities with respect to two 
norms (Lemma 11 and Corollary |12[), and an error estimate in a lower-order norm based on the 
Aubin-Nitsche trick (Lemma! 


Lemma 9. Let Assumption 1 hold true. Then, for z/ > 0, the discrete scheme Q is consistent. 
That is, the solution u of ^ satisfies 


Ak{u, v) = (/, v)r Vu G Xh. 


Proof. By Assumption 1, u € H'^(T) for an r > 1/2. In particular, u is continuous and vanishes 
on dr in the sense of traces. It follows that 


Ak{u,v) = (14curlru,curlru)r - /c^(14nu,nu)r + (Tfc(u), Vu G Xh. 

The integration-by-parts formula ([^ holds for (f := Vk curly-u and v G Xh (see pm for details 
concerning the Laplacian; they also apply to the Helmholtz case). The definition Q of Tk and 
relation Q then show the assertion. □ 

Lemma 10. There holds 


\Ak{v,w)\ <max{iy{s - 1/2) ^ (s - 1/2) ‘^}\\v\\H‘(r)\\M\H‘{r) yv,w € H^T), s > 1/2. 
Proof. By the continuity of Vk, and ( [l4) ) we obtain 

|(VfcCurlru,curlrrc)r| < || curlr u||^_i/2(r)ll curlric||j^-i/2(r) 


(s-1/2) 2|u|y^«(r)|'w|yfa(r) yv,w £ H%T), s > 1/2, 


and 


|( 14 nu,nu;)| < ||nu||f^_i 


jy-i/2(r)||nu;||^_i/2(r) < lkllL2(r)lhllL2(r) yv,w £ L‘^{T). (17) 


Combinations of (12) with the Cauchy-Schwarz and triangle inequalities, and estimate (13), 
yield 


KMi M) 7 l ~ (■® - 1/2) ^lbllH=(r)l|w^llH'>(r) yv,w £ s > 1/2 

and 

\{Tkiv), [u;]).y| -k \{[v],T_kiw))^\ <{s- l/2)~‘^\\v\\Hs(r)\\M\H‘{T) '^v,w £ s > 1/2. 

The previous bounds prove the assertion. □ 
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Lemma 11. Let Assumption 3 hold true. For any e > 0 there exists cq > 0 such that for 
V > h~^ there holds 


^ -cg(s- 1/2) '^h 


-2j,l-2s||^||2 




Proof. Application of ®, ©, and the fact that > 1 prove that there exists cg > 0 (G refers 
to Carding) such that 


|(Vocurl7-n,curl7-n)7— k‘^{Vknv,nv)r + i^{[v], [u]). 


> Ihlh 1/2 


CG||t'|li2(r) Vu G Xh. 


We are left with bounding the remaining terms. 


By Assumption 3, (14) and the inverse property we can bound 


(18) 


I ((H - M)) curlr V, curlr v)r\ < II curlr u||j^-i/ 2 (r) II curlr u||j^_i/ 2 -r,(r) 

<{s- 1/2 )~^\v\h‘>{t)\\ curlr u||j^_i/ 2 -.,(r) 

<{s- l/2)-'^h^/‘^-^\v\Hi/2(^-r)\\ curlr i;|lr-i/ 2 -rc(r), (19) 

for any v G and s G (1/2,1]. In order to estimate the last term above we use that there holds 


curlj Vj\\jj-i/ 2 -rc(j'^'i ^ Il'Cjllri/a-i-ccr.-u d ~ , J- 


40’ 


This follows from Fourier analysis, considering Fj as a sub-domain of and since ||'/’|lr‘(r ) ~ 
||(/°||r*(R 2 ) for t G [—1,1] and (j) G G“(Fj) with fp denoting ist extension by 0. By a standard 
domain decomposition estimate and bound © we then conclude that 


curlr u||^_i/ 2 - 


4r) 


< II curlr u||rV-l/2-r.^-^-^ 5, llt'l 


lr-i/ 2 -r'c(r) ~ l|c|iri/ 2 -rc(r) 


|u||j^l/2-rc (r) 


( 20 ) 


for any v G Xh- Combination of (19) and (20), and Young’s inequality, prove that 

|2 


( 21 ) 


|((H - M))curlru,curlru)r| </ 3 |^^|j:^i/ 2 (r )+^5 G (s-l/2) k iiv|i^i/2-rc(r) 

for any v G X^, s G (1/2,1], and p > 0. The two remaining terms are included analogously as 
in the case of the Laplacian [k = 0) considered in [U Lemma 4.4], Specifically, using (13) with 
s = 1/2 -|- e, Young’s inequality and the inverse property, one proves that 

\{TkV, [u]).^| -h |(H,r_fcu)^| < /i“^"^|u|^G2(r) + ^IINIIl 2 ( 7 ) '^veXh, 6>0, e> 0. (22) 


A combination of (18), (21) with p small enough, and (22) shows that there exist constants 
Cl, C 2 , cg > 0 (cg possibly different from before) such that 


\Mv,v)\> (l- cih ^'^)|v|^i/2(r) + y)llMlli2 (.^)-cg(s- 1/2) V 2"||u||^i/2-.,(r) 

for any v G Xh, 5 > 0, e > 0, s G (1/2,1], With <5 := selecting v > and replacing 4e by 
e we obtain the assertion. □ 
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Corollary 12. Let Assumption 3 hold true, and let e, e' > 0 be given. There exists cq > 0 such 
that for V > hT^ there holds 

\^k{v,v)\ > h^^'\\v\\]^y 2 (^r) ~ ^ ® ^ (1/2,1]- 

Proof. We select s = 1/2 + e' in (16) and use the inverse property to conclude that 

Il^^lli2(r) < /i"^'Vl^i/2(r) + z^llMlli2(^) Vu G Xh. 

This means that IIII^ 1 / 2 (- 7 -) ^ lkll^i/ 2 ^^^ for riny v G Xfi, and the assertion follows from 
Lemma [TTJ □ 


Lemma 13. Let Assumptions 1,2, and 3 hold true and assume that the meshes defining Xf, are 
globally quasi-uniform, i.e., h h. Given e > 0 choose v ~ /i“^. Then there exists hg > 0 such 
that the discrete scheme ([^ is uniquely solvable for h < ho. Furthermore, for h < ho there holds 

\\u - Uh\\H^-ri„o{r) ~ max{n{s - l/ 2 )"\ (s - l/ 2 )" 2 }||u - Uh\\H‘>{T) 

for any s G (1/2, min{riso, Treg}] - 

Proof. We first show the error estimate, i.e., for the time being let us assume that there is a 
(unique) solution Uh to Note that there holds Ar{v,w) = A-r{w,v) for any r G M and 
sufficiently smooth functions v, w. This follows from the fact the V-r is the adjoint operator 
of Vr. Let (j) G Lf^‘=°(r) be given (cf. Assumption 2). By standard approximation results there 
exists (fih G Xh such that 


\\4’ 4'h\\H‘‘{T) ~ ^ ll'/’lljTiso(r)’ 1/^ < -s < nso- (23) 

Using integration by parts (analogously to proving consistency in Lemma we hnd that there 
holds 


{u - Uh, lU_fc(/)r = {W-k4>, u-uh) = A_k{4>, u-Uh) = Ak{u - Uh, </>) = Ak{u -Uh,(j)- (fh)- 


Lemma 10 and (^23| then prove that 


|(u -'Uh,lT_fc(/>)r| < niax{z 2 (s - 1 / 2 ) b(s-l/ 2 ) ^}\\u - Uh\\H‘>(T)^ 


Ihhso (r)- 


Noting that || • || • by ( [IT| ) since r^so £ (1/2,1), this bound implies the 

error estimate via duality and by making use of Assumption 2 : 


~ sup 


|(k - Uh,'if)r\ 
IIV’llH’'iso-i(r) 


sup 

o/^effuso (r) 


Ku - Uh, W-k(f)r 
ll</ll ffuso (r) 


We are left with showing unique existence of Uh for small h. Since we are dealing with a quadratic 
discrete system, it is enough to show uniqueness. Therefore, in the remainder of this proof, we 
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assume that we are solving the homogeneous problem Q, i.e., / = 0 and u = 0. We have to 
show that only = 0 solves the homogeneous discrete scheme Q. The first part of this proof 
and the inverse property show that there holds 

ll4llHi-nso(r) < max{z/(s - l/2)-\ (s - l/2)-^}\\ul\\Hs(r) 

< /jri,o+i/2-2. - l/2)-\ (s - l/2)-^}\\ul^\\Hi/2^r) 

for s G (1/2,1], Now, if Hso < 1/2+ rc, then \\u^\\jji/ 2 -r^(j-^ < ||'W° ||j:/i-nso(r)> combination 
of the bound above with the estimate by Corollary |12| yields 

0 = > (/>=■' - P'l) 

for a number C{s) depending on s. For e,e' > 0 small enough, we can select s G (l/2, (hso + 
1 — e — e')/3) and find ho > 0 such that 

- C(s)/i2^‘-+2-6.-2. ^ ^2.' A _ ^^^^^2n,„+2-6s-2.-2.'\ ^ q ^ 


If ^’iso 

^1/2+rc 


> 1/2 + Tc then we additionally use the inverse property to bound 
Then we obtain, instead of (24), 



jyl/2-rc(7-) 


< 


(25) 


Analogously as before, for e, e' > 0 small enough, we can select s G (l/2,1/2 + (rc — e — e')/3) 
and hnd ho > 0 such that 


= h^^' (^1 - >0 Vfi < ho- 


In both cases, (24) respectively (25) proves that rtV = 0 for h sufficiently small. 


□ 


4.3 Proof of Theorem I 


The existence and uniqueness of Uh solving ([^ for small h is guaranteed by Lemma |13[ The 
proof of the error estimate follows the standard Strang strategy, that is, adding and subtracting 
a discrete function, using the triangle inequality, discrete Garding’s inequality (Corollary [l 2 |), 
consistency (Lemmaand boundedness (Lemma 10). More precisely, given e' > e > h small 
enough, v ~ /i“^, and v G X^, we hnd by the just mentioned arguments 


< ||u - -^11^1/2(7-) + (s - l/2)"^h^"^*"^"'||u/i - + h~‘^^'\Ak{u -v,Uh- u)| 

^ Iju — + (s — 1 / 2 ) h W'^h ~ v\\^l/2-rc(J-^ 

+ max{z 2 (s - l/ 2 )“\ (s - l/2)~‘^}{5~^\\u - v\\\s(j-) + (26) 
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for any <5 > 0. In the last step we applied Young’s and the inverse inequality. We have to 
consider the relation between riso and r^, cf. the proof of Lemma 13 In the case nso < 1/2 + Tc 
we bound \\uh - < \\uh - Otherwise, 




hr) 


\Wh- v\\jji/2-r,(^Y) - ''"°\\uh - v\ 

by the inverse property. Both cases are considered by 

\\uh - v\\jji/ 2 -r,^-j-) < h~°‘\\uh - 'y||Hi-nso(r)> “ ■= max{0,riso - Tc - 1/2}. 


We use the just established estimate in (26), and bound the norms of Uh — v by adding and 
subtracting u, applying the triangle inequality, and then bound — with the help 

of Lemma 13 This yields (with hidden constants depending on s) 

|2 




i//i/ 2 (r) 


_^^l-2s-2e -2au_ 


2 

H" ' ISO 


(T) 


+ h-‘^^'u{d-'\\u-v 


-ii 




+ 6h 


l-2s-2e' 


l^\\u fJ'hW jj-l/2(j-'^ 


5h" \\u — v\\ 

_|_ ^l+2riso-4s-2e'-2a,,2 


^ 11 '^ 


(27) 


The last term is handled yet again by the same technique (adding and subtracting v, inverse 
property of Uh — v): 


\u — 'W/i||^s(7-) < \\u — ^*(^11^4 — 'y||^i/2(-7-) + 11^^ — 'W/i||^i/2(-7-)^ • 


(28) 


Combination of (27) and (28), and reordering terms yields 

(^1 - ||m - Uh\\]ji/2^-j-j 


< /jl ^"||u — u||^l-r- (n-. 

~ II ''H ISO {y) 

,l-2s-2e\ 


+ (^1 + + h^+^nso-6s-2e’-2a^2^ 

+ u6~^ + /il+2nao-4s-2e'-2o^2^ 


u — V 


Hi/2(r) 


u — v\ 


H^iT) 


for two constants ci,C 2 > 0. We now select e' > e > 0 sufficiently small such that there is 
s G (l/2, (1 + riso — a — e — e')/3). This is possible since 1 + riso — a > 3/2. Furthermore, for 
the selected s we choose 6 = /i 2 ^f-i+e+ 2 e'+e" £q]. ^ Then the factor on the left-hand side is 

bounded from below by a positive constant for h being small enough. Replacing 6 and v also 
on the right-hand side shows that 


I Il2 ^ i,l-2s-2e'-2aii ||2 , ii ||2 

\U — Uh\\^\/2(J-^ \\U — + y-U — 


^^l-2s-(2e+4e'-He") ^l-H2ri,o-4s-2a-2e'-2e)^ | 


U — V\ 


HHT) 
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for h small enough. By the assumptions riso G (1/2,1) and G (0,1/2] one finds that, for e, e' 
sufficiently small, the term /;,i-2s-(2e+4e+e )||^ _ jg dominating one of the upper 

bound in the sense of best approximation orders in h. Then, renaming e + 2e' + ("/2 to be 
the new e/ this yields the error estimate of Theorem for the previously noted selection of 
s G (1/2, (1 + Tiso — a — e — e')/3). It is also clear that the upper bound for s can be dropped, 
as long as s < r^eg. 

5 Numerical results 

We consider the model problem Q with T = (—1/2,1/2) x (—1/2,1/2) x {0}, right-hand side 
function / = 1, and wave number k = 5. We use a decomposition of T into three sub-domains, 
as indicated in Fig.[^ and consider rectangular meshes which are piecewise uniform with respect 
to sub-domains, and globally quasi-uniform. The initial four meshes are also shown in Fig. 
The discrete spaces Xh consist of piecewise bilinear polynomials which are continuous on sub- 
domains. 


Figure 1: Decomposition of F into three sub-domains and initial mesh sequence. 

According to Corollary]^ and taking into account Remark]^ (i), we expect that for sufficiently 
large u, the error ||u — Uh\\jji/ 2 (^-i-'^ has convergence order close to 1/2, the optimal one for a 
conforming method and piecewise (bi)linear functions on quasi-uniform meshes, cf. [T]. However, 
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since the exact solution u of Q is unknown, the error cannot be computed directly. But even 
knowing u it would be difficult to calculate the necessary norm. For the Laplacian (k = 0) 
the residual and the L^-norm of the jumps form a reasonable upper bound for the error, see 
the discussion in m Section 5]. In the Helmholtz case (k ^ 0) energy arguments leading to 
such estimates do not apply without perturbation terms, see m Section 5]. Nevertheless, we 
conclude from the previously mentioned discussions that the sum of the two terms 

llullex - Re{{VkCurlrUh,c\irlrUh)T - k‘^{VkUh,Uh)r) 

'' -V-- 

= :residual =:jumps 

is a reasonably justified upper bound for the error ||u — Uh\\fji/2(^Y)- Here, ||u||gx is an approxi¬ 
mation of Re{WkU,u)r generated by extrapolation on a sequence of uniform meshes, cf. [8]. 

Figure shows the errors on a double logarithmic scale versus the inverse of the maximum 
over all side lengths. For comparison also the error in energy norm for the conforming variant on 
a sequence of uniform meshes and the curve 0.25/i^/^ are given. They confirm the convergence 
order of the conforming BEM. The results of the Nitsche approximation with = 

10 ,100,1000 indicate that, for v sufficiently large {u = 1000 in this case) this optimal order is 
achieved. At least for the model problem, this wave number and for the meshes considered, 
we do not observe a reduced convergence order. Such a reduced order can be seen in the case 
u = 10. For 12 = 100 the residual appears to reflect some pre-asymptotic behavior whereas the 
jumps still indicate a reduced convergence order. 

For illustration, we also present some conforming and Nitsche approximations to the solution 
u of Q, again with wave number k = 5. Figure]^ shows a conforming approximation (including 
homogeneous boundary condition) whereas Figure presents the Nitsche results for different 
meshes (the real parts on the left and the imaginary parts on the right). The coarser mesh 
(upper plots) is the last one from Figure 


l/Z 


+ ll[^/j]llL2(.y) 
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Figure 2: Errors of conforming BEM and Nitsche-BEM with v = 10,100,1000. 



Eigure 3: Conforming approximation for k = 5, real part (left) and imaginary part (right). 
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